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We report on a measurement of the T(1S+2S+3S) — > e + e~ cross section at midrapidity in p + p 
collisions at y/s = 200 GeV. We find the cross section to be 114 ± 38 (stat. + fit)^ 24 (syst.) pb. 
Perturbative QCD calculations at next-to-leading order in the Color Evaporation Model are in 
agreement with our measurement, while calculations in the Color Singlet Model underestimate it 
by 2a. Our result is consistent with the trend seen in world data as a function of the center-of-mass 
energy of the collision and extends the availability of T data to RHIC energies. The dielectron 
continuum in the invariant mass range near the T is also studied to obtain a combined cross section 
of Drell-Yan plus b-b — > e + e~ . 
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I. INTRODUCTION 

The main focus of the heavy flavor program at RHIC is 
to investigate the properties of the Quark-Gluon Plasma 
by studying its effect on open heavy flavor and quarkonia 
production. J / ip suppression induced by Debye screening 
of the static Quantum Chromo Dynamics (QCD) poten- 
tial between cc pairs was originally hailed as an unam- 
biguous signature of Quark-Gluon Plasma (QGP) forma- 
tion However, this simple picture is complicated by 
competing effects that either reduce the yield, such as 
co-mover absorption [2i J3j, or enhance it, such as in re- 
combination models |^-l6|- Recently, a growing interest 
in studying the T meson and its excited states has been 
kindled as it is expected that color screening will be the 
dominant effect contributing to any observed suppression 
of bottomonium production in heavy-ion collisions. A full 
spectroscopy of quarkonia states is now clearly recognized 
as one of the key measurements needed to understand the 
matter produced in high-energy heavy- ion collisions Q. 
In particular, it has been recognized that data on the 
particle spectra of bottomonia can provide valuable in- 
formation to constrain QGP models Q. Due to the low 
production cross section of b-b at RHIC (<r b _^ ral.9 /ib, 
[9J), recombination effects in A + A collisions are neg- 
ligible. At the same time, the interaction cross section 
of bottomonium with the abundantly produced hadrons 
in these collisions is small [lfj, so suppression due to 
absorption by hadronic co-movers is expected by these 
models to be relatively unimportant. However, it will 
still be important to study T production in d+Au colli- 
sions since available measurements by E772 [ll[ of cold 
nuclear matter effects on T production at lower energy 
show some suppression. Nevertheless, the amount of sup- 
pression seen for the T family is measured to be smaller 
than for charmonia. Therefore, bottomonium is expected 
to be a cleaner probe of high-temperature color screening 
effects. 

In addition to its important role in establishing decon- 
finement, a measurement of the T IS, 2S, and 3S states 
in p + p and heavy-ion collisions can help to set limits 
on the medium temperature. The quarkonium measure- 
ments help in reaching these key goals because (i) an 
observation of suppression of T production in heavy-ions 
relative to p + p would be a strong argument in support 
of Debye screening and therefore of deconfinement }12j |. 
and (m) the sequential suppression pattern of the ex- 
cited states is sensitive to the temperature reached in the 
medium 0- In this regard, lattice QCD studies have seen 
a burst of activity in recent years. Studies of quarkonia 
spectral functions and potential models based on lattice 



QCD indicate that while the T(3S) melts even before 
the deconfinement transition and the T(2S) is likely to 
melt at RHIC (Vs= 200 GeV), the T(1S) is expected 
to survive @, EH EH- Recent results [H, EH indicate 
further that almost all quarkonia states (J/tp, , Xc 
Xbi T(2S)) melt below 1.3 T c and the only one to sur- 
vive to higher temperature is the T(1S), which melts at 
2T C , where T c sw 175 MeV is the critical temperature 
for the parton-hadron phase transition. Therefore, a sys- 
tematic study of all quarkonia states in p + p, d+Au, and 
Au+Au collisions will provide a clearer understanding of 
the properties of the Quark-Gluon Plasma. 

Suppression of the T(2S) and T(3S) should be measur- 
able at RHIC energies with increased integrated luminos- 
ity. In the near future, the larger luminosities proposed 
by the RHIC II program [lj| will allow for a statistically 
significant measurement of all 3 states. With the ob- 
jective of embarking in such a long program, one of the 
first steps is to establish a baseline cross section measure- 
ment of the bottomonia states in p + p collisions. There 
are no previous measurements of T production in p + p 
at the top RHIC energy for heavy ions (an upper limit 
was estimated in the 2004 data with only half of the 
calorimeter |l7|). The luminosities available at RHIC 
in the 2006 run provided the first opportunity to mea- 
sure bottomonium at the previously unexplored center- 
of-mass energy of — 200 GeV. A dedicated trigger 
algorithm exploiting the capabilities of the STAR electro- 
magnetic calorimeter is essential for this measurement, 
and its development in STAR allows the T family to be 
studied in the e + e~ decay channel. In this paper, we 
report our result for the T(1S+2S+3S) cross section at 
midrapidity, obtained with the STAR detector in p + p 
collisions at y/s = 200 GeV via the e + e~ decay channel. 
This measurement uses an integrated luminosity of 7.9 
pb" 1 collected during RHIC Run VI (2006). We com- 
pare our data to perturbative QCD calculations done at 
next-to-leading order (NLO pQCD) in the Color Evapo- 
ration Model (CEM) [lj| and in the Color Singlet Model 
(CSM) pjj. 

The article is organized as follows. Section iHl explains 
the detector setup. The details of the quarkonia triggers 
are explained in Sec. Mil The acceptance and trigger 
efficiency are discussed in Sec. IIVI After a detailed dis- 
cussion of the data analysis procedure in Section [Vj we 
present our results and compare with pQCD calculations 
and with available data in Sec. lVIl Our findings are sum- 
marized in Sec. I VIII 



II. DETECTOR OVERVIEW 



* Deceased 



The main detectors used in the STAR quarkonia pro- 
gram are the Time Projection Chamber fTPO Hp . the 
Barrel Electromagnetic Calorimeter (BEMC) [2fJ, and 
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the Beam-Beam Counters (BBC) 

The BBCs are segmented scintillator rings covering the 
region 3.3 < \r)\ < 5.0. The STAR p + p minimum-bias 
trigger requires the coincidence of hits in the BBCs on 
opposite sides of the interaction region, and is used to 
monitor the integrated luminosity for measuring absolute 
cross sections. The BBC acceptance and efficiency were 
studied previously |22j | , and it was found from simulations 
that for p + p non-singly diffractive (NSD) events the 
BBC trigger was 87 ± 8% efficient. The absolute p + p 
cross section seen by the BBC trigger was measured via 
a van der Meer scan to be 26.1 ± 0.2 (stat.) ± 1.8 (syst.) 
mb. 

The TPC and BEMC have a large acceptance at midra- 
pidity: the single track coverage of the TPC+BEMC in 
pseudorapidity is |ry| < 1 and complete in azimuth. The 
BEMC is divided in r\ and (p into 4800 towers. The size 
of each tower in 77 x <p is 0.05 x 0.05 rad. The geomet- 
rical acceptance of the STAR BEMC for detecting both 
electrons from an T decay is about ~ 30% over all phase 
space, and « 60% for T's in the kinematic range \y\ < 0.5 
and p± < 5 GeV/c. We focus on this kinematic range, 
where the acceptance of the detector is optimal. 

The capability of the TPC+BEMC for electron iden- 
tification and the triggering capabilities of the BEMC 
are the two pillars of the STAR quarkonium program. In 
particular, the BEMC trigger allows us to sample the full 
luminosity delivered by RHIC to look for the high-mass 
dielectron signals characteristic of the T — > e + e~ decay. 



III. THE STAR QUARKONIA TRIGGERS 

The STAR quarkonia trigger is a two-stage system 
comprising a Level-0 (L0) hardware component [23| (de- 
cision time of w 1 /is) and a Level-2 (L2) software com- 
ponent (decision time of « 100 /js in p + p, and sa 400 fis 
in Au+Au). There were two separate T triggers used 
during the 2006 run; in the following we refer to them as 
Trigger I and Trigger II. They were identical in concept, 
with the only difference being an increase in thresholds, 
which was necessary to reduce the trigger rate and re- 
duce deadtime while data-taking. We discuss the triggers 
next, and list all the parameters and thresholds. These 
were chosen based on simulations and on expected trigger 
rates based on calorimeter data taken in 2003. 



A. T LO Trigger 

The L0 trigger is a minimum-bias p + p trigger with 
the additional requirement of signals in the BEMC con- 
sistent with a high energy electron. The energy deposited 
in the calorimeter towers is measured by collecting scin- 
tillation light from the electromagnetic shower. The sig- 
nal is digitized, pedestal subtracted, and sent to the L0 
Data Storage and Manipulation (DSM) hardware as a 6- 
bit ADC value for triggering. The calorimeter calibration 



is done such that tower ADC values are proportional to 
deposited electromagnetic transverse energy (E±), mak- 
ing it useful for triggering on high-pj^ electrons. The L0 
trigger decision is made for every bunch crossing. 

For the T analysis, the L0 decision was based on two 
quantities. The first is related to the signal in distinct 
towers. The tower with the highest E±, as represented by 
the DSM-ADC value, is called the "High Tower" (HT) 
in the event. The second quantity is a sum of towers 
in fixed rj-ip regions of the BEMC, where each region 
comprises 4x4 towers. Each of these regions is named a 
"Trigger Patch" (TP). Since each tower covers 0.05 units 
in 77 and 0.05 radians in 92, the TP coverage is 0.2 x 0.2 
rad. The L0 trigger is issued for Trigger I (Trigger II) 
if all of the following three conditions are met: (i) an 
event has a High Tower with a DSM-ADC value > 12 
(16), corresponding to E±m 2.6 (3.5) GeV deposited in 
the tower; (ii) the Trigger Patch containing the High 
Tower has a total DSM-ADC sum over the 16 towers in 
the patch with a value > 17 (19), corresponding to E± « 
3.8 (4.3) GeV, and (in) the STAR p + p minimum-bias 
trigger is met. The minimum-bias trigger is based on 
a BBC coincidence and is described elsewhere [22]. The 
coincidence of High- Tower and Trigger Patch triggers will 
be referred to as "HTTP" in the remainder of the paper. 

Figure [T] shows the p^ distribution for simulated elec- 
trons and positrons from T decays in which both daugh- 
ters are within the STAR acceptance. For each decay, we 
plot the p± distribution of the electron (filled circles) or 
positron (open squares) with the highest p± (the abscissa 
is simulated p±). The electron and positron distribu- 
tions are identical, as expected for a two-body decay into 
daughters with identical masses. The histogram in the 
figure is the sum, i.e., the p± distribution for the hard- 
est of the two daughters. The corresponding DSM-ADC 
values used in the L0 trigger (which are proportional to 
p±_) are shown in Fig. [5] The L0 HT Trigger II threshold 
of 16 is shown by the vertical dashed line. 

From the p± and ADC distributions for the hardest T 
daughter in Figs. [T] and [3J the average values are (p±) = 
5.6 GeV/c and (ADC) = 18.0 counts. Since the trigger 
threshold is placed at 16 ADC counts, the hardest daugh- 
ter will typically fire the trigger. We find that 25% of the 
T's produced at midrapidity have both daughters in the 
BEMC acceptance and at least one of them can fire the 
L0 trigger. The details of the HTTP trigger efficiency 
and acceptance are discussed in Sec. [IV] 



B. T L2 Trigger 

Once a L0 HTTP trigger is issued, the information for 
all detectors in STAR begins to be digitized. During this 
time, the L2 system can use the information from the 
4800 individual calorimeter towers to decide whether to 
keep the event or to abort the readout. The L2 system 
can use the full energy resolution of the calorimeter (a 
10-bit ADC value, in contrast to the 6-bit ADC used at 
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FIG. 1: The E± distribution for simulated T daughter elec- 
tron (filled circle) or positron (open square) with the highest 
E±. We show only daughters that fall in the STAR BEMC 
geometrical acceptance. The histogram is the sum of the two 
distributions. The L0 HT Trigger II threshold of 16 counts 
corresponds to E± ^3.5 GeV. 



FIG. 2: The L0 DSM-ADC distribution for the electron (filled 
circle) or positron (open square) with the highest ADC. The 
histogram is the sum. The simulated T is triggered when 
one of its daughter electrons is above the dashed line, which 
indicates the L0 HT Trigger II threshold of 16 counts. 



L0). The L2 T trigger makes the decision by looking 
for calorimeter signatures consistent with the production 
of a high invariant mass electron-positron pair. (Since 
the L2 trigger uses BEMC information only and does 
not distinguish electrons from positrons, we will refer to 
electrons and positrons simply as "electrons" in all dis- 
cussions of the L2 trigger.) At a minimum, one pair of 
candidate electrons is required. They are processed as 
follows. (We will denote kinematical quantities obtained 
in the L2 stage with the superscript "L2" to distinguish 
them from similar quantities obtained during the offline 
analysis.) 

The algorithm starts by searching for all towers above 
the L0 High Tower threshold. Each of these is treated as 
a seed for a 3-tower cluster. To produce the clusters, we 
search in the rj-ip region around each High Tower, with a 
search window of 3 x 3 towers. This area is smaller than 
the TP size in order to focus on electron finding, as our 
simulations show that electrons will likely have most of 
their energy contained in only 3 towers. We sort these 8 
surrounding towers according to their measured energy 
and pick the two high.est-.E_L towers in this list to pro- 
duce a 3-tower cluster. This L2 tower clustering gives 
a better estimate of the electron energy compared to a 
single tower due to possible shower leakage into neighbor- 
ing towers. If the energy of this 3-tower cluster is greater 
than 4.0 GeV it is considered for further processing, and 
we label such clusters as "L2 Cluster-1". 

We next look for additional electron candidates in the 



event. While it is possible for an event to have two tow- 
ers above the HT threshold, the majority of events will 
have only one. We select additional electron candidates 
by starting with towers which have E± > 0.7 GeV (10- 
bit ADC of at least 75 counts after pedestal subtraction). 
Starting from these second seeds, we again construct 3- 
tower clusters (L2 Cluster-2) using the procedure out- 
lined above. We require that the L2 Cluster-2 energy 
be E% 2 > 2.5 GeV. After the complete iteration to 
find electron clusters, we make all possible cluster pairs 
(where each pair must have at least one L2 Cluster-1) 
and calculate two pairwise quantities: the opening an- 
gle and the invariant mass of the pairs. These are cal- 
culated under the approximation that the vertex loca- 
tion is in the center of the detector and that the elec- 
trons travel in a straight line. Since the majority of the 
e + e~ pairs from T decays have a large opening angle 
(8), we look for pairs with cos(0f 2 2 ) < 0.5, where 6>f 2 2 
is the angle between the two clusters calculated in the 
L2 algorithm. With the combined information on the 
energy of the two clusters and their opening angle, we 
reconstruct the invariant mass via the approximate for- 
mula M[ 2 2 w y/2E^ 2 E^ 2 {\ - cos6»f 2 2 ). We can neglect 
the electron mass of 0.511 MeV/c 2 as it would only con- 
tribute s» 1 MeV/c 2 , which is small for ultra- relativistic 
electrons. For comparison, the straight-line approxima- 
tion and the energy resolution result in a mass resolution 
of « 900 MeV/c 2 as will be discussed in the following sec- 
tion. We select events with cluster-pair invariant masses 
in the range 6 < M^ 2 < 15 GeV/c 2 . If there is any pair 
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FIG. 3: The L0 DSM-ADC distribution (oc E±) for the high- 
est EMC tower of a candidate pair. We show Trigger I data 
after applying Trigger II thresholds (red circles) and Trigger 
II data (black squares). The yields are normalized by the 
integrated luminosity. The histograms are from simulation 
of T(1S), showing the corresponding distribution for electron 
daughters satisfying acceptance (solid line) and trigger (dot- 
ted line) requirements. The simulation histograms are nor- 
malized assuming B x da/dy— 100 pb, times a factor of 500 
for clarity. The vertical line is the Trigger II threshold of 16 
counts. The histograms have the bin centers set at integer 
values to match the integer nature of ADC counts. 



that satisfies both of the pairwise conditions in the event, 
the algorithm issues an L2 trigger. 



C. Trigger Performance 

In order to evaluate the trigger performance, we ap- 
plied the Trigger II cuts to the Trigger I data. To do 
this, the exact same trigger condition was applied offline 
on the recorded values of the original trigger input data. 
Since events satisfying Trigger II cuts also satisfy Trigger 
I cuts, applying the tighter set of Trigger II cuts to the 
Trigger I data offline allows us to make a single dataset 
with uniform properties for the entire 2006 run. The in- 
tegrated luminosity for Trigger I was C = 3.12 ± 0.22 
(syst.) pb -1 , and for Trigger II it was C = 4.76 ± 0.33 
(syst.) pb -1 , giving a total integrated luminosity of 7.9 
± 0.6 pb -1 , where the 7% uncertainty originates from the 
uncertainty in the BBC measurement of the cross section 
as determined by a van der Meer scan [22| . 

Figure [3] shows the ADC distribution of the tower with 
the largest E± for each candidate pair seen at the trigger 
level in an event (in very few cases we find more than 
one tower above the L0 threshold, so most events have 



only one candidate pair). The two triggered data sets are 
shown, with Trigger II displayed as squares and Trigger I 
(after applying Trigger II cuts) displayed as circles. The 
rejection factor achieved with Trigger II, defined as the 
number of minimum bias events counted by the trigger 
scalers divided by the number events where the trigger 
was issued, was found to be 1.8 x 10 5 . The distributions 
seen at the trigger level only include information from the 
BEMC. This causes the trigger rate seen in the experi- 
ment to be dominated by background from di-jet events 
with two nearly back-to-back 7r°s. This di-jet background 
is removed offline when including tracking information 
from the TPC. The trigger distributions are scaled by the 
overall luminosity in order to compare the relative nor- 
malization of the two datasets. The scale is chosen such 
that the counts in each triggered dataset correspond to an 
integrated luminosity of 1 pb -1 . The relative luminosity 
normalization between the datasets agrees to a level of 
~ 1%. The solid-line histogram is from simulated T(1S) 
after acceptance selection, requiring both electrons to de- 
posit at least 1 GeV of energy in the BEMC. It shows the 
spectrum for the daughter with the highest ADC count. 
The spectrum for the simulated T(1S) events that satisfy 
all the trigger requirements is shown as the dashed-line 
histogram. The vertical line at 16 ADC counts represents 
the L0 ADC threshold for Trigger II. In order to compare 
the size of the expected T signal relative to the trigger 
background from di-jcts, the histograms from the simula- 
tions for the T(1S) are scaled with two factors. The first 
factor corresponds to normalizing to an integrated lumi- 
nosity of 1 pb -1 assuming a cross section times branching 
ratio into e + e~ for the T(1S) at midrapidity of 100 pb, 
chosen because it is of the expected order of magnitude. 
From this, we expect one upsilon every w 3000 triggers, 
so we use a second multiplicative factor of 500 for dis- 
play purposes. The spectral shape of the data in Fig. [3] 
is therefore dominated by di-jet background. The same 
normalization and scale factors are used in Figs. HHSlfor 
comparing the trigger-level background distributions and 
the expected T(1S) signal. 

Figure [4] shows the distribution of the L2 Cluster- 1 
energy for all clusters that include a tower above the 
HT DSM-ADC threshold of 16 ADC and that have their 
corresponding Trigger Patch DSM-ADC sum above the 
threshold of 19 ADC counts. It can be seen that the peak 
of the L2-cluster distribution near 6 GeV is not right at 
the threshold of 4 GeV (vertical line) . This must happen 
because another trigger selection that is correlated to the 
measurement of the L2-cluster is being applied. In our 
case it is the L0 selection, which consists of both the 
HT and the TP requirement. Once we see a tower with 
energy above the HT threshold that is also in a trigger 
patch with energy above the TP threshold, the energy in 
the 3-tower cluster for that tower will likely already be 
above the L2 El threshold. 

The distribution for the L2 opening angle f5f 2 2 is shown 
in Fig. [SJ It is highly peaked towards back-to-back 
topologies, much more so than the distribution from sim- 
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FIG. 4: The distribution of L2 Cluster-1 energy Ei for all 
towers above the HT threshold in a trigger patch above the TP 
threshold. The L2 trigger requires Ei 2 > 4.0 GeV (vertical 
line). The line histograms show the T (from simulation) af- 
ter acceptance requirements (solid line), and after all trigger 
requirements (dashed line). The normalization and scaling 
factors are the same as in Fig. [3] 



ulated T's. This again reflects the fact that the triggered 
distribution is dominated by back-to-back 7r°'s. The ma- 
jority of these are rejected offline when requiring the pres- 
ence of a corresponding electron track in the TPC. 

The L2 invariant mass distribution shown in Fig. [5] is 
peaked at m 8 GcV/c 2 due to the cluster energy require- 
ments. The vertical lines depict the lower and upper 
thresholds at 6 and 15 GeV/c 2 . The histogram shows 
the simulated T(1S) before any trigger cuts (solid) and 
after passing all the previous thresholds (dashed). The 
trigger preferentially rejects lower L2-mass T events, be- 
cause these can happen when the energy clusters mea- 
sured in the calorimeter are lower than their average, 
and these clusters are preferentially rejected by the algo- 
rithm. From these simulations we estimate the T mass 
resolution of the L2 trigger to be 849 ± 8 MeV/c 2 . It is 
dominated by the BEMC energy resolution and by the 
straight-line approximation used to calculate the open- 
ing angle in the L2 trigger algorithm. This is about an 
order of magnitude larger than the offline mass resolu- 
tion, where the electron trajectories and momenta are 
obtained from tracking. Since all the T's that reach this 
stage are contained within the invariant-mass limits, this 
invariant mass cut serves mainly to reject background. 

In addition, the L2 information (such as that shown 
in Figs. |3HS]) is available online after every run. These 
distributions serve as useful diagnostic tools during data 
taking. 



FIG. 5: The distribution of cosine of L2 opening angle 9 12 for 
accepted events. The line histograms show the T distribution 
after acceptance requirements (solid line), and after all trig- 
ger requirements (dashed line). The vertical line shows the 
location of the trigger threshold. Normalization and scaling 
factors are the same as in Fig. [3] 



IV. T ACCEPTANCE AND TRIGGER 
EFFICIENCY 

To determine the geometrical acceptance of the detec- 
tor for measuring T — > e+e~, we combined two types 
of events in the following way. For the first type, we 
performed GEANT simulations of the T — > e + e~ de- 
cays and propagated them through the detector geome- 
try. The simulations were done with uniform population 
of the p±-y phase space, folded with a Gaussian in y 
(a = 1) and a realistic p± distribution. We chose the 
form dN/dp± oc p±/ (exp(p±/T) + 1) with the parameter 
T = 2.2 GeV/c obtained from a fit to CDF data 0| . The 
dependence of the acceptance on the choice of p± distri- 
bution is negligible. This was verified by using functional 
forms derived from data at lower energy. Each simulated 
decay was combined with a simulated p + p minimum- 
bias event using the PYTHIA event generator [25[ with 
CDF Tune A settings (26|. For the second type, we used 
a set of data collected by triggering on a RHIC bunch 
crossing. These were labeled "zero-bias" as they do not 
require signals in any of the STAR detectors. These zero- 
bias events do not always have a collision in the given 
bunch crossing, but given that in the 2006 p + p run the 
dominant contribution to the TPC occupancy was from 
pileup, these events will include the pileup from out-of- 
time collisions and all detector effects. When the zero- 
bias events are combined with simulated events, they pro- 
vide the most realistic environment to study the detec- 
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FIG. 6: The L2 invariant mass M 12 distribution for accepted 
events. The line histograms show the T M^-i distribution 
after acceptance requirements (solid line), and after all trig- 
ger requirements (L0 HTTP trigger, L2 cluster, L2 opening 
angle and L2 invariant mass) of the e + e~ daughters (dashed 
line). Vertical lines show location of the trigger thresholds. 
Normalization and scaling factors are the same as in Fig. [3] 



FIG. 7: Combined acceptance x efficiency as a function of ra- 
pidity yr for the decay T— > e + e~ in STAR. Down triangles: 
BEMC acceptance only; open squares: adding L0 require- 
ment; upright triangles: adding L2 requirement; open circles: 
adding TPC acceptance+TPC tracking cuts and TPC-BEMC 
matching. 



tor efficiency and acceptance. Each PYTHIA+T event 
was embedded into a zero-bias event. The vertex posi- 
tion chosen for the simulated event was sampled from a 
realistic distribution of event vertex positions obtained 
from the T triggered data. With this procedure, we es- 
timate both the trigger efficiency and the reconstruction 
efficiency (discussed in Sec. IV) . 

For our purposes, a simulated T is considered to fall 
in the acceptance of the detector if each of its decay elec- 
trons deposits at least 1 GeV in a BEMC tower, as sim- 
ulated by GEANT. We find that the acceptance depends 
strongly on the T rapidity (j/x), but weakly on its trans- 
verse momentum (p±x)- At |yx| < 0.5 the acceptance 
is 57%, dropping to below 30% beyond \yy\ — 0.5, and 
is close to zero beyond |yx| = 1-0. This is illustrated 
in Fig. [71 where the downward triangles depict the geo- 
metrical acceptance of the BEMC detector e geo obtained 
from the analysis of the simulated T decays in the real 
data events. 

For those simulated T's that are accepted, we calcu- 
late the trigger efficiencies. This requires simulations of 
the BEMC response, digitization, and running the result 
through the offline reconstruction software chain. The 
HTTP requirement and the L2 trigger condition cut out 
an additional fraction of the events, and the combined 
acceptance and trigger efficiency is shown in Fig. [7J The 
squares include the effect of the HTTP trigger require- 
ment £lo and the upright triangles include the L2 trigger 



TABLE I: Acceptance, trigger efficiency, and tracking effi- 
ciency for reconstructing T — > e + e~ in STAR in the kinematic 
region \yr \ < 0.5. The first 4 rows are for the IS and the last 
row is for the cross section-weighted 1S+2S+3S combination. 
All systematic uncertainties are listed in Table [TV] 



Quantity Value 

e geo 0.57 

Cgco x clo 0.25 

Cgco X CL0 X 6L2 0.21 

Cgco X 6L0 X £ L 2 X Ctrack X £R 0.13 

Egco X 6L0 X £L2 X e t rack X E R (1S + 2S + 3S) 0.14 



condition £l2- The open circles include TPC acceptance 
plus tracking efficiency etrack and TPC-BEMC matching 
6r. These selections, as well as the particle identification 
cuts, are discussed in Sec. |V] For the region |yx| < 0.5, 
the acceptance and efficiencies are given in Table HI 

The relative systematic uncertainty on the L0 HT trig- 
ger efficiency was estimated to be and is the domi- 
nant source of systematic uncertainty of the T trigger effi- 
ciency. The systematic uncertainties for the L0 efficiency 
are asymmetric because the underlying E±_ distribution 
of the daughter electrons is not fiat in the threshold re- 
gion (see Fig. [3]). A shift to a higher threshold reduces the 
yield more than a shift to a lower threshold increases it. 
The L0 HTTP trigger has largest effect on the efficiency, 
as only 43% of the events in the acceptance remain after 
the trigger condition is applied. After including the L2 
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TABLE II: 


Branching fractions for T(nS) — > 


e+e" [23 and 


total cross 


sections at y/s = 200 GeV from 


an NLO CEM 


model [JJ|. 


See text for details. 




T state 


B (%) 


a (nb) 


T(1S) 


2.38±0.11 


6.60 


T(2S) 


1.91±0.16 


2.18 


T(3S) 


2.18±0.21 


1.32 



trigger, tracking, and TPC-BEMC matching efficiencies, 
we obtain a combined efficiency of 13.2% for the T(1S). 

A similar procedure was applied to the T(2S) and 
T(3S) states, which have slightly higher efficiencies due 
to their larger masses. We calculated a weighted aver- 
age among the 3 states, including the branching ratio 
and the ratio of cross sections, in order to obtain an 
average efficiency to be applied to the measured yield 
of T(1S+2S+3S). For this we use the branching ratios 
compiled by the PDG 



27] 



and cross section ratios from 
NLO pQCD calculations in the Color Evaporation Model 
(CEM). The calculation used a bottom quark mass of 
m b = 4.75 GeV/c 2 , the PDFs used are MRST HO [Hj], 
the choice of scale is \i = tut, and the center-of-mass 



energy is y/s = 200 GeV |16|. The branching ratios and 
cross sections we used for this purpose are shown in Ta- 
ble HI 

This procedure results in a combined acceptance, 
trigger efficiency, track finding efficiency, and TPC- 
BEMC matching efficiency of 14.3% for the averaged 
T(1S+2S+3S) combination, see the last row of Table Q] 
We estimate the sensitivity of the efficiency to the values 
in Table ILTl in two ways, (i) By varying the branching ra- 
tio values by their uncertainty, we find that the efficiency 
varies in the range 14.2-14.4%. (ii) Instead of using the 
numbers from Table ILTl we can use measured values of the 
ratios T(2S)/T(1S) and T(3S)/T(1S) from Ref. 0. No 
further input is needed because the weighted average of 
the efficiency does not depend on the overall scale of the 
cross section: it cancels out when averaging. By varying 
these measured ratios within their uncertainty we obtain 
weighted efficiencies in the range 14.1-14.4%. Therefore, 
our results do not depend strongly on the values in Table 
ILTl In particular, the overall scale of the cross section 
from the NLO calculations has no direct impact on our 
results. Additional analysis cuts related to electron iden- 
tification are discussed in the following section. 

With thresholds of « 4 GeV, the hadron rejection 
power of the BEMC towers at the energy relevant for T 
decay daughters is e/h~ 100. The main source of back- 
ground for the trigger comes from high-pj^ 7r°'s decaying 
into two photons that deposit energy in the BEMC. The 
BEMC will therefore trigger on two high-pj^ 7r°'s with a 
large opening angle. These events are typically from di- 
jets. Since these are also produced by large-momentum- 
transfer events with low cross section, the trigger rate 
is sufficiently small even in the presence of this back- 
ground, that the T trigger can sample all collision events 



and is limited only by luminosity. The triggered distribu- 
tions shown in Figs. [5HH1 are dominated by the ir° photon 
background. These background events are rejected in the 
offline analysis. 



V. T ANALYSIS 

During the offline analysis, we use a complete emula- 
tion of the trigger to obtain all candidate BEMC tower 
pairs for an event. We use the TPC to select charged 
tracks and require that they point close in rj-ip to the posi- 
tion of the candidate clusters. The TPC also allows us to 
obtain improved electron kinematics compared to those 
derived from BEMC information available at the trigger 
level and to perform particle identification via specific 
ionization. The matching of TPC tracks to BEMC clus- 
ters is also useful for electron identification via the ratio 
E/p. 

Tracks are selected based on the number of TPC points 
found during the track pattern recognition and used in 
the fit to obtain the track kinematic parameters. The 
TPC tracks can have a maximum of up to 45 space 
points. We require that all tracks have at least 52% 
of their maximum possible points, and a minimum of 
20 points. The first requirement guarantees we have no 
split tracks, and the second requirement sets a floor for 
those tracks which have a number of points smaller than 
45 due to passing through inactive areas of the detector. 
We select tracks with p± > 0.2 GeV/c to reject most 
low momentum tracks in a first pass. We do not im- 
pose higher p± requirements in the track selection stage 
because we also require E/p matching with the BEMC, 
where E > 2.5 GeV, as discussed below. Tracks are also 
required to extrapolate back to a primary vertex found in 
the event. We find the combined TPC acceptance times 
tracking efficiency for detecting each T daughter to be 
£tpc = 85%, and to be approximately independent of 
electron p±. Hence, the pairwise efficiency can be ob- 
tained by squaring the single particle efficiency. 

In order to guarantee that the analysis only uses tracks 
that could have fired the trigger, we impose a requirement 
that the tracks extrapolate close to the BEMC candi- 
date towers. The requirement used was for them to be 
within a circle of radius R < 0.04 in rj-<p space of the 
L2 candidate clusters. We find that this cut has an ef- 
ficiency of = 93% for a given T daughter electron, 
and has good background rejection power. This is il- 
lustrated in Fig. [HI which shows the R distribution for 
simulated T(1S) daughters as the line histogram. The R 
distribution from tracks passing basic quality cuts in the 
triggered events is shown as the solid circles. 

An important component of the analysis is the vertex- 
finding efficiency. We find that in contrast to minimum- 
bias p + p events, the T triggered events have a very high 
vertex-finding efficiency. The analysis benefits from the 
presence of the HT trigger, since these events are likely 
to have a high-p^ track that facilitates the task of find- 
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FIG. 8: The radial distance R = \/ At] 2 + Aip 2 between the 
TPC track and the EMC electron cluster for simulated T 
daughters (line histogram, arbitrary scale), and the R distri- 
bution from data (circles). We reject candidate tracks with 
R > 0.04 (vertical line). 



ing the vertex. We find that the vertex-finding efficiency 
£vcrtcx for T events is 96.4 ±0.9%. 

In the 2006 run, the luminosity was high enough that 
there can be multiple primary vertices due to pileup 
events in the TPC. We find that about w 9% of the 
T-triggered events have 2 or more vertices. For this 
analysis, we searched for candidates from all the vertices 
found in an event. We chose the vertex by requiring 
that the vertex also matched the high momentum TPC 
tracks that were already selected based on the BEMC 
tower extrapolation. Since the BEMC is read out after 
every bunch crossing, out-of-bunch pileup (interactions 
that happen before or after the triggered bunch cross- 
ing) is rejected by the TPC track-BEMC cluster match- 
ing requirement. We end up with a unique, unambiguous 
vertex in all events. Within-bunch pileup (multiple inter- 
actions in the same bunch crossing) is negligible for this 
analysis. Pileup rejection will become more important as 
the luminosity of RHIC is increased. 

Electrons were identified by selecting charged parti- 
cle tracks with specific ionization energy loss dE/dx in 
the TPC consistent with the ionization expected for elec- 
trons. In the momentum region of interest (p > 3 
GeV/c), there is overlap between electrons and charged 
pions, so a cut was placed that yielded an electron ef- 
ficiency of edE/dx = 84%. The cut was chosen to opti- 
mize the effective signal S e g of single electrons, S e ff = 
S/(2B/S +1), where S is the electron signal and B 
is the hadron background. To do this, we construct 
a normal distribution of ionization measurements. We 
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FIG. 9: E/p distribution from T-triggered data events, select- 
ing a high-purity sample of positrons via dE/dx (filled circles). 
We compare this to simulations from reconstructed e + and e~ 
from T decays (solid line) and from ty° events where the de- 
cay photons produce a conversion e + e~ pair (dashed line). 
The scaled sum of these two contributions is shown as the 
dot-dashed line. 



use the measured ionization of each track, compare it 
to the expected ionization for an electron and divide by 
the expected dE/dx resolution (which depends on the 
track length and number of measurements). This yields a 
normalized Gaussian distribution of ionization measure- 
ments, no dE/dx- We fit this distribution with one Gaus- 
sian function to represent the electrons signal S, and two 
Gaussian functions to model the background B from pi- 
ons and other hadrons. 

In addition, we used the combined information of 
the EMC energy E, as obtained by the 3-tower can- 
didate clusters from L2, and the TPC track momen- 
tum p to compute the E/p ratio, which should be unity 
for ultra-relativistic electrons to a very good approxi- 
mation. Figure IHl shows the E/p distribution compari- 
son from an electron sample obtained from the Trigger 
I and II datasets (filled circles). Cuts were placed at 
E/p = 1.0 ± 0.3, which is close to a 2er cut, given that 
fitting with a Gaussian gives a = 0.157 ± 0.002. A shape 
that is approximately Gaussian is expected for the E/p 
distribution due to the Gaussian shape of the resolution 
in both the energy measurement done in the BEMC and 
in the curvature measurement done in the TPC (curva- 
ture is proportional to l/p±). We see that the E/p distri- 
bution in the data shows a non-Gaussian tail. We studied 
the shape in several ways. To rule out distortions from 
misidentification of hadrons, we show only positrons to 
avoid any antiproton contamination and use tight dE/dx 
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cuts to select only a very high-purity electron sample. 
We studied the energy dependence (to rule out thresh- 
old effects near 4 and 2.5 GeV), isolation cuts on the 
electron tracks (to rule out contamination of the BEMC 
energy measurement from nearby charged particles), and 
the E/p shape of embedded electrons and photons. We 
conclude that the tail at high E/p has two sources: 

1. Electron bremsstrahlung. The solid line in Fig. [S] 
shows the simulated E/p for T electrons combined 
with a PYTHIA underlying p + p event, embedded 
into zero-bias-triggered events. Many of these pro- 
duce a Gaussian E/p distribution, but we also see 
a tail. By selecting simulated electrons which un- 
dergo bremsstrahlung, we find cases in which the 
Monte Carlo momentum at the outermost TPC 
point differs from the momentum at the T decay 
vertex by more than 100 MeV/c. We then look at 
the E/p of these electrons, and find that their mean 
is shifted. All entries in the region E/p ss 1.3 of the 
solid line in Fig. |H] come from these cases. There- 
fore, part of the non-Gaussian shape seen in the 
data is accounted for by electron bremsstrahlung 
in the detector material, and is included in the dis- 
tribution of simulated T daughters. 

2. Photon conversions. Events with a high-Q 2 di-jet 
that include back-to-back 7r°'s also fire the T trig- 
ger. Some of the tt° daughter photons will con- 
vert into e + e~ pairs and leave a TPC track point- 
ing at the EMC cluster. Generally, the e + and 
e~ will strike the calorimeter near the sibling pho- 
ton, resulting in a track with a high E/p value. 
We analyzed these simulated 7r° events and applied 
the same tracking cuts, calorimeter clustering, and 
BEMC-TPC matching used in this paper to the 
resulting electrons. The E/p distribution for the 
electrons in these simulated events is shown as the 
dashed line in Fig. |9l and we see a non-Gaussian 
tail extending to larger values than that of the T 
electrons. The average E/p for the electrons from 
7T° events is « 1.18 compared to ~ 1.08 for T elec- 
trons. 

Of these two effects, the first needs to be taken into ac- 
count in the T efficiency because the T daughters which 
undergo bremsstrahlung and have a resulting E/p value 
outside our cut will be removed from the analysis. The 
second effect does not need to be included in the E/p ef- 
ficiency since it is not due to bottomonium events. These 
7r°-7-conversion events will appear in the invariant mass 
distributions, but they can be subtracted as follows. A 
single photon conversion produces an unlike-sign pair. 
When there are multiple photon conversions in an event, 
the combinations that can be made include unlike-sign 
and like-sign electron pairs. The unlike-sign pairs from 
a real photon will have zero invariant mass and do not 
affect the analysis. The additional unlike-sign combina- 
torial pairs will have a distribution that can be modeled 



by the like-sign combinations and are therefore removed 
when subtracting the like-sign pair combinatorial back- 
ground. 

Since the E/p distribution in the T-triggered dataset 
has both types of events, we reproduce the shape of the 
data using the two types of electron simulations: the 
electron E/p distributions from T events (including the 
bremsstrahlung tail) and from ir° events. These are 
added and scaled to approximately fit the data in the 
region 0.6 < E/p < 1.7, shown as the dot-dashed line in 
Fig. M 

As discussed above, we need to determine the efficiency 
for electrons only from T sources. To do this, we use 
the T embedding simulations to estimate the efficiency 
of our E/p cut, which includes both the Gaussian shape 
expected from detector resolution and the non-Gaussian 
tail due to bremsstrahlung. The systematic uncertainty 
in the determination of the E/p cut efficiency is domi- 
nated by the uncertainty in the knowledge of the detector 
material (and hence on the amount of bremsstrahlung), 
which was estimated to be a factor of 2. The material un- 
certainty is different from that quoted in previous STAR 
electron analyses because in this paper we do not restrict 
the event vertex to be near the center of the detector, 
where the material budget is lower and is better known. 
To estimate the uncertainty in the E/p cut efficiency, we 
construct an E/p distribution by splitting the original 
embedding simulation sample into electrons that undergo 
bremsstrahlung and electrons that do not. We scale the 
number of those that do by a factor of 2 and sum the two 
distributions to obtain a new E/p distribution. The esti- 
mated efficiency found in the bremsstrahlung-augmented 
distribution is lower by sa 3%, which we assign as the 
systematic uncertainty of e B / p . All the efficiencies used 
to obtain the T cross section and their systematic uncer- 
tainties are collected in Table HVl below. 



VI. RESULTS AND DISCUSSION 

Once we have tracks that satisfy our criteria and are 
matched to tower clusters that satisfy all the L0 and 
L2 trigger requirements, we form electron-positron pairs 
to produce the invariant mass spectrum, using electron 
track momenta reconstructed in the TPC. We use the 
like-sign combinations of e + e + and e~e~ to estimate the 
combinatorial background (see e.g. [3(|) via 

= 2^N ++ N^ ■ A+ ~ (1) 
V A ++ A -- 

where denotes the unlike-sign pair differential in- 
variant mass distribution dN^ /dm and denotes 

the acceptance for unlike-sign pairs (similarly for the 
like-sign distributions). The symmetry of the BEMC 
and TPC for accepting unlike-sign and like-sign pairs 
makes the ratio of acceptances unity, so we only use 
2yJ N ++ N-- for the combinatorial background. The 
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FIG. 10: Unlike-sign raw yield N+- (filled circles) in the re- 
gion \y e &\ < 0.5, where y ee is the pairwise rapidity, and like- 
sign combinatorial background 2^/N ++ N (line histogram). 

The background is fit to an exponential multiplied with an er- 
ror function representing the turn-on of the STAR T trigger 
(dashed line). 



unlike-sign and like-sign background invariant mass spec- 
tra are shown in Fig. [TU1 



A. T Line Shape 

In this analysis, we cannot resolve individual states 
of the T family due to the limited statistics, finite mo- 
mentum resolution, and electron bremsstrahlung result- 
ing from the large material budget during the 2006 run. 
Therefore, the yield reported here is for the combined 
T(1S+2S+3S) states. The total dielectron unlike-sign 
yield after like-sign background subtraction has three 
contributions. 

1. The T states. 

2. The Drell-Yan continuum. 

3. The b-b continuum. 

In order to separate these three contributions, we ob- 
tain parametrizations for the expected shapes of each of 
the three components of the dielectron unlike-sign yield. 
The contribution from the T states is obtained by the 
same simulations that were used to obtain the detec- 
tor efficiency. The reconstructed invariant mass shape 
for the IS, 2S, and 3S states are individually obtained. 
Each of the shapes from the simulation is fit with a 
functional form introduced by the Crystal-Ball experi- 
ment [3 1], which can accommodate detector resolution 



FIG. 11: The e + e~ signal after subtracting the like-sign com- 
binatorial background. The solid line is the function used 
to fit the data, composed of: (a) three Crystal-Ball functions 
(Eq.[2]) that represent the combined T(1S+2S+3S) line shape, 
(b) the Drell-Yan, and (c) the b-b contributions. We fit the 
data using the integral of the fit function in each mass bin, 
shown by the dashed-line histogram. The sum of the two con- 
tinuum contributions is shown by the dot-dashed line. The 
integral of the Crystal-Ball functions provides the net T yield. 



and losses due to bremsstrahlung in the detector mate- 
rial, and has the form: 



f(m) = N 



exp(- 
A- (B- m=H) 



(m-fj.) \ 
m—fi\— n 



for 
for 



> -a 
< -a 



(2) 



Requiring that the function and its derivative be contin- 
uous constrains the constants A and B to be: 



.4 : 
B 



a 



(3) 



We fix the parameters of the three Crystal-Ball functions 
representing the three T states and then adjust the rela- 
tive scales according to the average branching ratios and 
according to the ratios of cross sections shown in Ta- 
ble HIl The overall integral of the Crystal-Ball functions 
combined this way is left as a free parameter in the fit, 
which allows us to obtain the T yield. 

Figure [IT] shows the data used to obtain the T and 
continuum yields. The data points are the unlike-sign 
e + e~ signal after subtracting the like-sign combinatorial 
background. The background subtraction is done bin- by- 
bin. The error bars shown are statistical. The fit includes 
the contributions from the T(1S+2S+3S) states and the 
continuum contributions from Drell-Yan and b-b shown 
as the solid-line function. The T states are modeled with 
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Eq.[2]to account for detector resolution and for losses due 
to bremsstrahlung which result in a tail to low values of 
invariant mass. The fit is done using the integral of the 
fit function in each bin, and is shown as the dashcd-linc 
histogram. We discuss next the procedure used to extract 
the continuum contributions (dot-dashed line) and the 
T(1S+2S+3S) yield. 



B. Drell-Yan and b-b Continuum Contributions 

The Drell-Yan continuum is parameterized from a 
pQCD NLO calculation [32| done in the kinematic range 
\yr\ < 0.5 and m > 4 GeV/c 2 . We convolute the cal- 
culated spectrum with the detector resolution (account- 
ing for bremsstrahlung), but this introduces very small 
changes in the shape of the spectrum. We find that the 
shape of the Drell-Yan continuum is well described by a 
function of the form 



A 



(1 + m/m ) r 



(4) 



with the parameters mo = 2.70 GeV/c 2 and n = 4.59. 

The b-b contribution is p arameterized from a simula- 
tion using PYTHIA 8 [25|. We turn on production of 
b-b pairs, and follow their fragmentation and subsequent 
decays of the B hadrons to look for dielectrons origi- 
nating from the b quarks. We convolute the simulated 
shape with the detector resolution and bremsstrahlung, 
and find that the shape is well described by a function of 
the form 



Am b 



(1 + m/m ) c 



(5) 



29.7 GeV/c 2 , and 



with the parameters b — 1.59, tuq 
c = 26.6. 

Since the STAR T trigger is meant to reject events with 
low invariant mass, we parametrize the trigger response 
on both continuum contributions by multiplying them 
with an error function 



erf((m — m tr i g )/w) + 1 



(6) 



where m tr i g is related to the trigger thresholds and w 
describes the width of the turn-on of the error function 
due to finite detector resolution. We obtain the parame- 
ters from the like-sign data (Fig. IT0|) by multiplying the 
error function with an exponential exp(— m/T) to ac- 
count for the random like-sign combinations at higher 
mass. We use the like-sign data since the trigger turn-on 
shape for like-sign and unlike-sign pairs is the same in 
our detector, but the like-sign invariant mass distribu- 
tion is purely due to combinatorics and can be fit with 
the simple parametrization given above (Fig. 1101 dashed 
line). The parameters are mtrig = 8.1 ± 0.8 GeV/c 2 and 
w = 1.8 ±0.5 GeV/c 2 . 
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FIG. 12: The 1-er (horizontal lines) and 2-a (diagonal lines) 
contours obtained by fitting the data, showing the favored 
regions for the Drell-Yan and b-b cross sections. The dashed 
line shows the prediction for the Drell-Yan cross section from 
a pQCD NLO calculation j32j. 



To obtain the T yields and cross section, we perform 
a fit of the unlike-sign invariant mass distribution after 
subtracting the like-sign background including contribu- 
tions for the T(lS-r-2S+3S), and the continuum due to 
Drell-Yan and b-b. Since the extracted T yield from the 
fit will be sensitive to the continuum yield, we next dis- 
cuss the effect that variations of the Drell-Yan and b-b 
cross sections can have on our result. 

We find that the resulting shapes of the two continuum 
contributions are very similar. We can fit the continuum 
yield using the parameterized Drell-Yan and b-b contri- 
butions. The 1- and 2-er contours in the 2-D parameter 
space <J b _i vs. ctdy are shown in Fig. [T21 where the cross 
sections quoted are corrected for efficiency and accep- 
tance, and measured in the phase space region \y\ < 0.5 
and 8 < m ee < 11 GeV/c 2 . 

The fit shows a strong anti-correlation between the 
Drell-Yan and the b-b cross sections due to the simi- 
larity of the parameterized shapes. The shape of the 
Drell-Yan is slightly favored by our data, and we ob- 
tain the minimal \ 2 f° r a Drell-Yan cross section of ss 
38 pb with a negligible b-b contribution. The fit gives 
X 2 /NDF = 1.1. The prediction from an NLO calculation 
of the Drell-Yan cross section [32[ is shown as the ver- 
tical line at 38.6 pb, and is consistent with values in the 
1-er region (horizontal lines). Our data also allow values 
of the b-b cross section up to s» 15 pb within the 1-er 
range, provided the Drell-Yan yield is reduced. Given 
the anti-correlation, our data are mainly sensitive to the 
sum of the two cross sections, and within the 1-er con- 
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TABLE III: Extraction of T(1S+2S+3S) yield by bin count- 
ing and fitting. The sums are done in the range 8 < m < 11 
GeV/c 2 . 







iV+_ -2s/N ++ N— 


80.9 


T counts 


59 ±20 


Bin-by-bin Counting 




N+- - 2ViV ++ JV__ 


82.7 


T counts 


61 ± 20 


Single-bin Counting 




N+- - 2y/N++N— 


75 ±20 


T counts 


54 ± 18 



tour the allowed range for the sum of cross sections is 
(cdy ± 0&-&)||y|<o.5, 8<m<n GcV/c 2 = 38 ± 24 pb. This 
variation is taken into account in the fit to extract the T 
yield, and is included in the quoted uncertainty. 



C. T yield and cross section 

The fitting function in Fig. ITT1 includes a total of three 
free parameters: the T yield, the Drell-Yan yield and 
the b-b yield. This allows us to extract the T yield di- 
rectly, automatically including (i) the statistical uncer- 
tainty for each mass bin in Fig. II H and (it) the uncer- 
tainty due to the anti-correlation between the continuum 
contributions and the T contribution. To compare the 
yield extraction using a different method, one can also 
perform bin counting. There are two ways to do this. 
One method is to take the background-subtracted unlikc- 
sign yield directly from the data shown in Fig. [TTJ where 
the background subtraction is done for each mass bin. 
We then sum the resulting histogram in the mass re- 
gion of the T peak, 8-11 GeV/c 2 . This sum includes 
the T yield and the continuum contribution, so we sub- 
tract from it the contribution from the Drell-Yan and b-b 
continuum obtained in the fit. A second method is to 

sum the yield of the unlike-sign (A r _ ] ) in the region 8- 

11 GeV/c 2 , do the same for the like-sign positive (N ++ ), 
and like-sign negative (N ), and then do the subtrac- 
tion N + - — 2y/N ++ N In other words, in the first 

method we do the subtraction bin-by-bin and then do 
the sum, in the second method we do the sum first to get 
a single bin and then we do the subtraction. The results 

for estimating — N ++ N-- from the combined 

fit, the bin-by-bin counting method, and the single-bin 
counting method, are shown in Table IIIII 

The quoted uncertainty in — 2 ^N ++ N is listed 

for the single-bin counting method, which is obtained 
from a straightforward application of the statistical er- 
rors of the corresponding yields (A^_| = 230, N ++ = 92, 

A^__ = 65). It amounts to a contribution of 26% to the 
uncertainty, and originates only from counting statistics. 
In addition to this purely statistical uncertainty, we must 
also take into account the uncertainty of the continuum 



subtraction. This part of the uncertainty is obtained 
from the fitting method in all the quoted T counts in 
the Table. It should be noted that one cannot use count- 
ing statistics to estimate the uncertainty on the T yield 
due to the continuum subtraction. The reason is that 
the T and continuum yields are obtained from the same 
data, namely the invariant mass dielectron spectrum in 
Fig. [TT] The two yields are therefore anti-correlated: a 
larger continuum yield reduces the extracted T counts. 
One must take this anti-correlation into account in the es- 
timation of the uncertainty on the T yield. The best way 
to do this is to use the fitting method. In the fit, both the 
statistical precision of the invariant mass spectrum and 
the anti-correlation between the T and continuum yields 
in the estimation of the uncertainty are taken into ac- 
count when varying the T yield and the continuum yield 
to minimize the \ 2 - We find the uncertainty from the fit 
to be 33% of the T yield, and we hence use this fitting- 
method result to estimate the uncertainty in all the T 
yields quoted in Table [TTTT We quote this as the "stat. ± 
fit" uncertainty in our cross section result. It should be 
noted that with the statistics of the present analysis, we 
find that the allowed range of variation of the continuum 
yield in the fit is still dominated by the statistical error 
bars of the invariant mass distribution, and so the size 
of the 33% uncertainty is mainly statistical in nature. 
However, we prefer to denote the uncertainty as "stat. 
± fit" to clarify that it includes the estimate of the anti- 
correlation between the T and continuum yields obtained 
by the fitting method. A systematic uncertainty due to 
the continuum subtraction can be estimated by varying 
the model used to produce the continuum contribution 
from b-b. These variations produce a negligible change in 
the extracted yield with the current statistics. 

When comparing different counting methods in Ta- 
ble IIIII we see that the difference between the bin-by- 
bin counting method and the fitting method is negligible. 
The single-bin counting method yield is lower than the 
one from the bin- by-bin counting method by 9%, and we 
assign this as a component of the systematic uncertainty. 

Table GID also lists the yield of T(1S±2S±3S). In the 
fitting method, the yield and the quoted uncertainty is 
obtained directly from the fit, as well as the contribu- 
tion from Drell-Yan and b-b discussed in the previous sec- 
tion. This has the advantage that the uncertainty due to 
the continuum subtraction and the correlations between 
the T, Drell-Yan and b-b contributions are automatically 
taken into account when exploring the parameter space. 
The disadvantage is that there is a model dependence on 
the line-shapes used for the fit. In the counting methods, 
to extract the T yield we subtract the continuum con- 
tribution obtained from the fit. This reduces the model 
dependence on the T line shape. However, we must still 
account for the uncertainty in the estimate of the con- 
tinuum contribution in the determination of the T yield 
uncertainty. Since this uncertainty should include simi- 
lar correlations between the T yield and continuum yields 
as found in the fitting method, the relative uncertainties 
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should be approximately equal. Therefore, we use the 
same relative T yield uncertainty for the counting meth- 
ods as for the fitting method. To get the total T yield we 
must correct the above numbers for the yield outside the 
integration region. This correction can be obtained from 
the fitted Crystal-Ball functions, and gives an additional 
9% contribution to the T yield. We report results for the 
cross section using the bin-by-bin counting method. 

In order to transform the measured yield of 
T(1S+2S+3S) into a cross section, we applied several 
correction factors: 



#(nS) x cr(nS) 



N 



n=l 



Ay x e x C ' 



(7) 



where the symbols are as follows. B(nS) is the branch- 
ing fraction for T(nS) — > e + e~. a(nS) is the cross sec- 
tion da/dy for the nS state in the region \yx\ < 0.5. 
N = 67 ± 22 (stat.) is the measured T(1S+2S+3S) 
yield from the bin-by-bin counting method in Table IIIII 
with a 9% correction to account for the yield outside 
8 < m ee < 11 GeV/c 2 . Ay = 1.0 is the rapidity interval 
for our kinematic region \yr\ < 0.5. The total efficiency 
for reconstructing members of the T family is the prod- 
uct £ = e gco X Cyertex X £L0 X £L2 X e T PC X £r X f-dE/dx x e E / p , 

where the symbols are as follows. e gco is the BEMC ge- 
ometrical acceptance. e ver tex is the vertex-finding effi- 
ciency. £lo and cl2 are the trigger efficiencies for L0 and 
L2, respectively, ctpc is the TPC geometrical acceptance 
times tracking efficiency for reconstructing both daugh- 
ters in the TPC. e R is the TPC-BEMC rj-tp matching 
efficiency. £dE/dx is the electron identification efficiency 
from the specific ionization requirement, and £e/p is the 
electron identification efficiency from the E/p selection. 

We find for the cross section at midrapidity in y/s = 
200 GeV p + p collisions the result 



3 

E 

71=1 



B(nS) x cr(nS) 



114 ±38 tf 4 P b 



(8) 



The uncertainties quoted are the 33% statistical±fit un- 
certainty (mentioned in the discussion of Table IIII[) and 
the systematic uncertainty, respectively. 

The major contributions to the systematic uncer- 
tainty are: the uncertainty in the choice of bin-counting 
method, the uncertainty in the integrated luminosity, the 
uncertainty in the BBC efficiency for p + p NSD events 
and the uncertainty in the L0 trigger efficiency for T 
events. The polarization of the T states also affects the 
estimation of the geometrical acceptance. We estimate 
this uncertainty by comparing simulations of fully lon- 
gitudinal and fully transverse decays and comparing the 
acceptance of these cases with the unpolarized case. A 
list of all corrections and systematic uncertainties in the 
procedure to extract the cross section is compiled in Ta- 
ble [IVJ The combined systematic uncertainty is obtained 
by adding all the sources in quadrature. Note that the 



TABLE IV: Systematic uncertainties on the measurement of 
the T cross section. 



Quantity 



Value Syst. uncertainty on da/dy (%) 



iV+_ -2^N++N— 


82.7 


+0 
-9 


C 


7.9 pb" 1 


±7 


EBBC 


0.87 


±9 


£geo 


0.57 


+3.0 
-1.7 


^vertex 


0.96 


±1.0 


£L0 


0.43 


+7.5 

-5.9 


£L2 


0.85 


+0.7 
-0.2 


ETPC 


0.85 2 


2 X ±5.8 


C-R 


0.93 2 


2X + 1 ' 1 


(■dE/dx 


0.84 2 


2 x ±2.4 


e E / P 


0.93 2 


2 x ±3.0 


Combined 




+ 22 8 Db 

-24.1 P u 



single-particle efficiencies enter quadratically when re- 
constructing dielectron pairs, so we multiply the single- 
particle uncertainty by a factor of 2 when estimating the 
pairwise uncertainty. 
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FIG. 13: The STAR measurement of the midrapidity 
T(1S±2S±3S) cross section times branching ratio into elec- 
trons (star). Error bars are statistical, the box shows the 
systematic uncertainty, and the scale is given by the left axis. 
The raw yield vs. y is shown by the histogram at the bottom 
(diagonal-line fill pattern), with scale on the right axis. The 
cross section was calculated from the yield between the verti- 
cal dot-dashed lines, \y e e\ < 0.5. The open squares are from 
an NLO CEM calculation, and the two dotted lines give the 
limits for the prediction from a NLO CSM calculation of the 
T cross section (see text). 

The result we obtain for the cross section is shown in 
Fig. 1131 where the datum point given by the star symbol 
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is our measurement, the error bars and the box depict 
the statistical and systematic uncertainties, respectively. 
To illustrate the acceptance in rapidity, we also show the 
unlike-sign pairs after like-sign background subtraction, 

AT + _ - 2^/N ++ N , in the T region 8 < m ee < 11 

GeV/c 2 as a hashed histogram. The scale on the right 
axis of the figure is used for the counts in the histogram, 
and the scale in the left axis of the figure is used for the 
cross section. We compare our measurement with NLO 
CEM predictions [32| of the T(1S) rapidity distribution. 
Since we measure all three states and only in the dielec- 
tron channel, the calculation of the T(1S) is scaled by a 
factor 



B(1S) x ct(1S) + B(2S) x cr(2S) + B(3S) x ct(3S) 



(9) 



in order to compare it to our measurement of the cross 
section for all three states. The branching ratios and 
cross sections used for this scale factors are those from 
Table HU The calculation is in agreement with our mea- 
surement. The two dotted lines in the plot are the upper 
and lower bounds of the cross section obtained from a 
calculation in the CSM for direct T(1S) production [l8[ 
based on NLO code developed for quarkonium produc- 
tion at hadron colliders 33]. Since the calculation is for 
the IS state alone and for direct T production (ignoring 
feed-down from P-states), to compare to our measure- 
ment, which includes all 3 states and feed-down contri- 
butions, the values from the calculation were divided by 
a factor 0.42 to account for this (see Ref [l8[ for details). 
The bounds in the calculation are obtained by varying 
the bottom quark mass and the renormalization and fac- 
torization scales. The CSM prediction is lower than our 
data, indicating that additional contributions are needed 
beside production via color singlet. 

In Fig. [TH we also compare our T(1S+2S+3S) result 
with measurements done in p + A, p + p and p + p col- 
lisions at center-of-mass energies ranging from 20 GeV 
up to 1.8 TeV [H, HI, IH-[4l| , and to NLO CEM predic- 



tions 16] for a wide range of center-of-mass energies. 

Our result is consistent with the overall trend, and 
provides a reference for bottomonium production at the 
top RHIC energy. 



VII. CONCLUSIONS 

The STAR experiment has measured the 
T(1S+2S+3S) — > e + e~ cross section at midrapid- 
ity, \y ee \ < 0.5, in p + p collisions at ^fs — 200 GeV to be 
(B x da/dy) ls + 2S + 3S = 114 ± 38 (stat. + fit) +24 (syst.) 
pb. Calculations done in the Color Evaporation Model 
at NLO are in agreement with our measurement, while 
calculations in the Color Singlet Model underestimate 
our cross section by « 2a. Our result is consistent 
with the trend as a function of center-of-mass energy 
based on data from other experiments. We report 
a combined continuum cross section, Drell-Yan plus 
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FIG. 14: Evolution of the T (1S+2S+3S) cross section with 
center-of-mass energy for the world data and an NLO CEM 
calculation. The error bars on the STAR datum point are 
statistical and systematic as in Fig. 1131 



b-b — > e + e , measured in the kinematic range \y e e\ < 0.5 
and 8 < m ee < 11 GeV/c 2 , of (cr DY + a b _ h ) = 38 ± 24 pb. 
The STAR measurement presented here will be used as a 
baseline for studying cold and hot nuclear matter effects 
in d+Au and Au+Au collisions, as the relatively clean 
environment provided by the STAR high-mass dielectron 
trigger permits the approach outlined in this paper to 
be deployed up to the most central Au+Au collisions. 
With increased luminosity, a better determination of 
the cross section, its p± dependence and a separation 
of the 2S and 3S states will be possible. The projected 
luminosity upgrades to RHIC should increase the T yield 
to rj 8300 in p + p and ~ 11200 in Au+Au collisions 
during one RHIC year [l6|. The increased statistics will 
greatly reduce the uncertainty in the determination of 
the continuum cross section and will allow a thorough 
study of the bottomonium sector by resolving the 2S 
and 3S states. 
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